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Abstract 

We analyze the breakdown of Fermi-liquid behavior within the 2D Hubbard 
model as function of doping using our recently developed numerical method 
for the self consistent summation of bubble and ladder diagrams. For larger 
doping concentrations the system behaves like a conventional Fermi-liquid 
and for intermediate doping similar to a marginal Fermi-liquid. However, for 
smaller doping pronounced deviations from both pictures occur which are due 
to the increasing importance of the short range antiferromagnetic spin fluc- 
tuations. This is closely related to the experimental observed shadow states 
in the normal state of high- Tc superconductors. Furthermore, we discuss the 
implications of our results for transport experiments. 
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Central questions in the theory of the high-Tc superconductors are related to the Fermi- 
liquid (FL) properties. Although the existence of a Fermi surface (FS) is well establishedi, 
pronounced deviations from FL behavior could be observed. For strongly overdoped com- 
pounds the experimental results seem to be consistent with a standard FL picture, whereas 
for systems near the optimal doping this is clearly not the case: indications for a linear 
energy-dependence of the inverse quasi particle lifetime reported in angle resolved photoe- 
mission (ARPES) measurement s0, the unusual temperature dependence of transport coef- 
ficientslii and the uj/T scaling behavior of the dynamical magnetic susceptibility observed 
in neutron scattering measurement^! can not be understood within the conventional FL 
picture. For underdoped systems the deviations from FL behavior are even stronger due to 
the dominance of short range antiferromagnetic correlations. To achieve a comprehensive 
understanding of the normal state properties of the high-Tc superconductors the marginal 
Fermi-liquid (MFL) theory was proposedll. In this phenomeno logical approach the charac- 
teristic energy scale is given only by the temperature and the inverse lifetime of the quasi 
particles is proportional to their excitation energy. However, this ansatz does not take into 
account the strong short range antiferromagnetic correlations in the CUO2 planes, whose 
influence on the quasi particle properties was demonstrated in NMrI, neutron scattering 
and also in ARPES measurementJli. 

Interesting insight into the structure of these highly correlated electron systems could be 
obtained from Monte Carlo and exact diagonalization studied!. An alternative approach to 
investigate the dynamical properties of the high-T^ systems is offered by perturbation theory. 
To achieve an appropriate description of the low energy spin fluctuations we use in the fol- 
lowing the fluctuation exchange approximation (FLEX) by Bickers and Scalapino0lll. This 
self consistent summation of all bubble and ladder diagrams is a conserving approximation 
in the sense of Baym and KadanofH. In particular, it takes into account the interaction 
of the quasi particles with spin fluctuations which are expected to be the dominating low 
energy excitations in the high-Tc materials. In distinction to other diagrammatic methods 
the self consistency allows to consider larger interaction strengths, whereas in distinction to 



exact diagonalization calculations and quantum Monte Carlo simulations we are able to use 
much larger lattices which gives access to the low energy properties with negligible finite 
size effects. 

Since we are interested in dynamical properties (especially in the quasi particle lifetimes), 
we transform the FLEX-equations, which are formulated on the imaginary frequency axis, to 
the real axis0. Thus we determine by calculating the momentum dependent Greens function 
G'k(co') and its self-energy Ek(to') the low energy properties within the framework of the 2D 
Hubbard model. This enables us to study the character of the elementary excitations as 
a function of the doping concentration x = 1 — n, where n is the occupation number per 
site. We present new results on (i) the continuous breakdown of the Fermi-liquid picture 
for decreasing doping (ii) the interplay between this phenomenon and the occurrence of 
shadows of the Fermi-surface and (iii) the temperature scale set by the spin fluctuations and 
its influence on the temperature dependence of the Hall-coefficient. 

We consider the one-band Hubbard Hamiltonian 



where c]^ is the creation operator of an electron at lattice site i and with spin a. t^j is 
the hopping matrix element between sites i and j, /i is the chemical potential and U is 
the local Coulomb repulsion. The results presented below are obtained for nearest neigh- 
bor hopping t = 0.25 eV and U = 4t. The self-energy of this Hamiltonian within the 
FLEX-approximation, neglecting the particle-particle excitations which were shown to be of 
minor importanceEl, is given by the following momentum and Matsubara-frequency su m§ 
T,k(iujn) = jj Z]k',n' Gk' ii^n')Vk-k/ {iujn~iujn') where T is the temperature and the number 
of sites of the finite lattice. The effective interaction Vk(^^'m) resulting from the summation 
of bubble and electron-hole ladder diagrams is given by 
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Here, Xq(^'^m) = -jjE: 
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'n 



3 



{2n + l)7rT and z/^ = 2m7rT are the fermionic and bosonic Matsubara frequencies, re- 
spectively. Furthermore, the Greens function is given by the Dyson equation Gk(iu;„) = 
{iun + fJ' — £^(k) — Sk(«ci;„))~^ where £:(k) = —2t{cos{kx) + cos{ky)) is the free dispersion. 
These equations are analytically transformed to the real frequency axis yielding a set of 
equations for the Greens function G^^uj) and the self-energy Sk(i-^), which is solved self 
consistently^. 

In Fig.l we present our results for a system with doping x = 0.10 at T = 63K. Fig. 1(a) 
shows ImEk(cc') for momenta close to the FS (solid circles in the schematic plot of the 
Brillouin zone in the inset of Fig.lc). Note that uj = corresponds to the Fermi energy. For 
low energy hole {u < 0) as well as for low energy particle {u > 0) excitations a linear in- 
dependence for energies \uj\ > uo' where cj' ^ 8 meV can be seen. Note that uo' is of the order 
of the temperature broadening T = 63i^ = 6 meV. Although this behavior is consistent 
with the MFL theory the different slopes for hole excitations and for particle excitations are 
in distinction to MFL. In Fig. 1(c) we plot the corresponding real part of the self-energy and 
in Fig. 1(e) the spectral density, which consists only of a main peak. 

In Fig. 1(b) we plot ImSk(co') for momenta k close to the Fermi surface shadow (FSS), 
which are indicated as open circles in the inset of Fig. 1(c). The FSS consists of those points 
k' of the Brillouin zone which are shifted by Q = (tt, tt) relative to the FS. In a recent papeJil 
we discussed the formation of shadow states resulting from short range antiferromagnetic 
correlations. Due to the strong coupling of quasi particle states with momenta k + Q at 
the FSS to those near k at the main FS spectral weight is transferred from k to k + Q 
which results in the formation of shadow states. We found this coupling to be mostly 
pronounced for FS and FSS states nearest to the (tt, 0)-point, resulting from the flat bands 
in this region. At the FSS, the deformations of ImSk(u;) relative to MFL behavior are 
much stronger compared to FS states and a double peak structure develops, which is mostly 
pronounced for k = (7r,7r/8) (indicated as (II) in the inset of Fig. 1(c)). In Fig.l(d) we plot 
the real part of the self-energy for the same momenta on the FSS. The local maxima close to 
the FS are clearly visible, but not strong enough to produce new poles of the Greens function, 
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which would result in new antiferromagnetically induced quasi particle states as proposed 
by Kampf and SchriefferS. Nevertheless, satellites in the density of states occur very close 
to the Fermi energy (Fig.l(f)), which build up the experimentally observed shadows of the 
Fsi. For momenta far away from the FS as well as from its shadow we find the FL-like 
quadratic frequency dependence of ImSk(co') (not shown). From the k-dependence of the 
spin susceptibility we calculate a rough estimation of the correlation length ^ by neglecting 
vertex correctionJll. For x = 0.10 we find ^ ~ 3.5ao, where denotes the lattice spacing in 
the Cu02 planes, in agreement with the experimental observations. 

For larger doping, the correlation effects in the high- Tc compounds are expected to 
get weaker resulting in conventional FL behavior. Therefore, we investigated the doping 
dependence of the self-energy. In Fig. 2 we show lmT,i^{u) for a FS state and for a FSS state 
for X = 0.10 and x = 0.20. For the FS momentum (solid lines) in both cases there exists a 
crossover from a quadratic to a linear tu-dependence. However, for x = 0.10 the crossover 
takes place at uj ^ 8meV (which is of the order of the temperature, see discussion above), 
whereas for x = 0.20 it happens at the much higher value ~ 25 meV. For increasing doping 
this crossover takes place at even higher energies resulting in FL behavior for x > 0.25 (not 
shown). Now we discuss the frequency dependence at the FSS (dashed lines). For x = 0.20 
the spectral density pk{uj) is small for u close to the Fermi level because shadow states arise 
only for low doping concentrations. Consequently, states shifted by Q relative to the FS 
exhibit the frequency dependence ImSk(co') oc u"^ like all other states which are far away 
from the FS as can be seen in Fig.2(b). However, for x = 0.10 (Fig.2(a)) strong deviations 
from this behavior occur due to antiferromagnetic fluctuations, as discussed in the preceding 
paragraph. 

The new structures of the self-energy, especially the double well in Im Sk(co') for momenta 
close to the FSS, show that in contrast to the MFL theory the temperature is not the 
only energy scale in this model. A second characteristic energy scale is set by the spin 
fluctuations themselves. To investigate the interplay of these two energy scales in more 
detail we solved our equations for a fixed doping x = 0.11 and various temperatures. In 
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Fig. 3(a) we show ImSk(Lij) for a momentum close to the FS (sohd circle labeled (I) in 
the inset of Fig. 1(c)). For increasing temperature the linear cu-dependence is replaced by 
a functional form like lmTiy_{uj) oc (u; — c^o)^ in an energy range of the order of T. The 
different slopes of the inverse quasi particle lifetime for electron respectively hole excitations 
are for higher temperatures reflected in a shift Uq of the maximum of ImSk(c(j) away from 
the Fermi energy. Furthermore, an extrapolation to even lower temperatures indicates that 
ImSk(c<J = 0) as function of T does not saturate to for T — > as in the FL picture 
but to a finite value. Such a behavior is consistent with our picture that the shadow states 
arise from quasi particle decay due to the strong coupling of states near to the FS to states 
near to its shadow. Of course, this coupling is also present at T = 0. In Fig. 3(b) we show 
ImSk(c<j) for the same doping and for a momentum close to the FSS (open circle labeled 
(II) in the inset of Fig. 1(c)). The double well structure vanishes at T ^ 300K. Note that 
ImSk(cL' = 0) does not only deviate from conventional FL or MFL behavior, but even has a 
maximum for T ^ 300K and afterwards decreases for decreasing temperature, again yielding 
for T — > a finite value. For states far away from the FS and from its shadow we detected 
a quadratic temperature dependence (not shown). The existence of a new energy scale can 
also be seen in the momentum averaged density of states (not shown). Here, the spin density 
wave pseudogap, which is clearly visible for low temperatures (T = 75K), becomes smaller 
with increasing temperatures and vanishes at T ^ 750K. Therefore, we conclude that this 
temperature represents for x = 0.11 the characteristic energy of the spin fiuct national. 

The unusual momentum, frequency and temperature dependence of the elementary ex- 
citations has important consequences for the transport properties of the high-Tc supercon- 
ductors. The linear tu- dependence of ImSk(co') for the FS states (see discussion of Fig.l(a)) 
corresponds to a linear temperature dependence of the electrical resistivity^ for materi- 
als near the optimal doping concentration in agreement with the experiments. However, 
for smaller doping we found deviations from the linearity of ImH^luj). Moreover, shadow 
states with very unconventional properties of the scattering rate get occupied that can ad- 
ditionally influence the transport coefficients. This results in deviations from the linear 



temperature behavior of the electrical resistivity. Furthermore, due to the coexistence of 
different temperature dependences for different scattering modes the Hall coefficient be- 
comes T-dependenti'0. To demonstrate that our self consistent FLEX-calculation is able 
to describe such an effect we plot in the inset of Fig. 3(b) our results for the Hall coefficient 
obtained by neglecting further vertex corrections for the conductivity within the magnetic 
fieldii for doping x = 0.11 and various temperatures. Rh{T) decreases for increasing tem- 
perature and saturates above the characteristic spin fluctuation energy of T 750K in 
agreement with the experimentJi. 

In conclusion, we investigated the properties of the elementary excitations in the normal 
state of high-Tc superconductors in the framework of FLEX and the one-band Hubbard 
model. Due to our new numerical method we calculated the self-energy directly on the real 
frequency axis and obtained new results for the flne structure of this quantity as function of 
doping. We found its frequency and temperature behavior to be strongly momentum depen- 
dent and to exhibit very unconventional structures especially for states at the Fermi surface 
shadow, which get more pronounced for lower doping. We obtained a continuous crossover 
from Fermi-liquid behavior for very large doping (x > 0.25) to an intermediate range which 
behaves similar to a marginal Fermi-liquid (x ~ 0.15) and flnally to a doping range where 
the system is dominated by strong antiferromagnetic correlations (x ~ 0.1). For increasing 
temperature we flnd in all doping regions a transition to Fermi-liquid behavior. The doping 
and temperature dependence of the Fermi-liquid properties is schematically illustrated in 
Fig. 4. Note, that the phase diagram of our theory is in agreement with the normal state 
phase diagram which is obtained by analyzing the transport properties of {La, Sr)2Cu04^E 
In a forthcoming paperB we investigate the consequences of this qualitatively new behavior 
of the quasi particle states for superconductivity. 
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FIGURES 

FIG. 1. Self-energy and spectral density for x = 0.10 and T = 63if. (a) ImSk(^^) for four 
different momenta which are close to the Fermi-surface. Note the linearity up to very small fre- 
quencies, (b) The same quantity for momenta close to the Fermi surface shadow. The double well 
structure is mostly pronounced for k nearest to the (tt, 0)-point. (c) ReSk('^) for momenta close 
to the Fermi surface, (d) The same quantity for momenta close to the Fermi surface shadow, (e) 
Spectral density along the main Fermi surface, (f) Spectral density along the Fermi surface shadow. 
Inset of (c) Brillouin zone for orientation. Momenta on the main Fermi surface are indicated by 
solid circles, those on the shadow Fermi surface by open circles. 

FIG. 2. Imaginary part of the self-energy for fixed temperature T = 63K and two doping 
concentrations, (a) ImSk(w) for small doping (x = 0.10). The solid line indicates a momentum on 
the main Fermi surface (solid circle labeled (I) in the inset of Fig. 1(c)), the dashed line a momentum 
on its shadow (open circle labeled (II) in the inset of Fig. 1(c)). Note the strong deformations relative 
to Fermi-liquid and marginal Fermi-liquid behavior at the shadow, (b) ImEk(w) for larger doping 
(x = 0.20) with solid and dashed line as in (a). The frequency dependence of both is Fermi-liquid 
like. 

FIG. 3. Imaginary part of the self-energy or a; = 0.11 and different temperatures. In (a) we 
plotted ImSk(a;) for a momentum close to the Fermi surface, in (b) for a momentum close to the 
Fermi surface shadow. The double well structure vanishes T k. 300i^. Note that the scattering 
rates are large for both cases and for all temperatures we considered. In the inset of (b) we show 
the Hall-coefficient for the same doping as a function of T. 
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FIG. 4. Schematic diagram for the transition from Fermi-hquid to non Fermi-hquid behavior 
in the normal states for intermediate couphng U = 4t as discussed in the text. The boundaries of 
the various areas are diffuse and there is a continuous crossover between the three regions. The 
smahest doping for which our calculations were done is x = 0.09. Therefore, the antiferromagnetic 
phase transition observed in the cuprates for even lower doping is not included in the diagram. 
Note that our results are consistent with the Mermin- Wagner theoremil because we found no 
antiferromagnetic instability. 
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